Correlations between vibrational entropy and dynamics in liquids 
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A relation between vibrational entropy and particles mean square displacement is derived in super- 
cooled liquids, assuming that the main effect of temperature changes is to rescale the vibrational 
spectrum. Deviations from this relation, in particular due to the presence of a Boson peak whose 
shape and frequency changes with temperature, are estimated. Using observations of the short-time 
dynamics in liquids of various fragility, it is argued that (i) if the crystal entropy is significantly 
smaller than the liquid entropy at T g , the extrapolation of the vibrational entropy leads to the 
correlation Tk ~ To, where Tk is the Kauzmann temperature and To is the temperature extracted 
from the Vogel-Fulcher fit of the viscosity, (ii) The jump in specific heat associated with vibrational 
entropy is very small for strong liquids, and increases with fragility. The analysis suggests that these 
correlations stem from the stiffening of the Boson peak under cooling, underlying the importance of 
this phenomenon on the dynamical arrest. 

PACS numbers: 64.70.Pf,65.20.+w.77.22.-d 



When a liquid is cooled sufficiently rapidly to avoid 
crystallization, the relaxation time r below which it be- 
haves as a solid increases up to the glass transition tem- 
perature T g where the liquid falls out of equilibrium. In 
strong liquids r displays an Arrhenius dependence on 
temperature, but in other liquids, said to be fragile, the 
slowing-down of the dynamics is much more pronounced. 
In general the dynamics is well captured by the Vogel- 
Fulcher law log(r) = C + U/(T - T Q ), although non- 
diverging functional forms can also reproduce the dynam- 
ics well As the temperature evolves, two quantities 
appear to be good predictors of r: the space available 
for the rattling of the particles on the picosecond time 
scales 0, Q , embodied by the particles mean square dis- 
placement (u 2 ) observable in scattering experiments, and 
the difference between the liquid and the crystal entropy 
When extrapolated below T g this quantity vanishes 
at some Tk, the Kauzmann temperature which ap- 
pears to correspond rather well to the temperature T 
extracted from the Vogel-Fulcher law Q . The correlation 
between dynamics and thermodynamics has been inter- 
preted early on as the signature of a thermodynamical 
transition at Tk toward an ideal glass where the coiffig- 
urational entropy associated with the number of meta- 
stable states visited by the dynamics, or inherent struc- 
tures, would vanish J3- This view is appealing and still 
influential today [a, p , although it is not devoid of con- 
ceptual problems [lof^ Elastic models fi~lT - fl3j propose an 
alternative scenario of the glass transition: fragile liquids 
are simply those which stiffen under cooling, reducing the 
particle mean square displacement and increasing the ac- 
tivation barriers that must be overcome to flow. In this 
view the rapid change of entropy in fragile liquids stem 
from the temperature dependence of the high-frequency 
elastic moduli [141 ] . This is consistent with some obser- 
vations supporting that fragile liquids stiffen more under 
cooling [12j ] , but predicts entropy variations in the liquid 
several times stronger than those observed [l3| . 



Such distinct interpretations of the liquids entropy still 
coexist because the latter consists of several contributions 
hard to disentangle experimentally [lj| , in particular the 
configurational and the vibrational entropies. Phrased 
in the context of energy landscape, the vibrational en- 
tropy corresponds to the volume of phase space associ- 
ated with a typical configuration 4j . The corresponding 
specific heat differs from the one of the phonon bath of a 
harmonic elastic network, both due to the non-linearity 
of the interactions, and to the fact that the vibrational 
spectrum of the inherent structure evolves with tempera- 
ture. The latter effect does not appear in high frequency 
dynamic heat capacities studies, as it requires changes of 
configurations that take place on the time scale r to oc- 
cur, but does contribute to the jump of specific heat char- 
acterizing the glass transition, whose amplitude is known 
to strongly correlate to the liquid fragility 16]. The con- 
figurational fraction of the jump in specific heat has been 
estimated using calorimetry in quenched and annealed 
glasses [l5j . measurements of elastic moduli [I7j . densi- 
ties of states (1814201 ] and non-linear dielectric suscepti- 
bilities [21]. Overall, a great variation was found among 
glasses, with a configurational fraction ranging from 15% 
to 80% and apparently decaying with the liquid fragility 
[2l| . Nevertheless, data are yet sparse and more analysis 
needs to be done to quantify respectively the correlations 
between dynamics and the different contributions to the 
entropy. In this Letter the relation between vibrational 
entropy and particle mean square displacement is inves- 
tigated. From available data on the latter quantity, it is 
argued that the vibrational entropy by itself displays the 
evoked correlations with the dynamics, and that these 
correlations must be induced by the stiffening of some 
soft degrees of freedom under cooling, rather than via an 
overall rescaling of the vibrational spectrum. 

On time scale t << r one may approximate a super- 
cooled liquid as a solid in a well-defined configuration. 
Scattering experiments yield a spectral analysis of the 
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corresponding dynamics, enabling to define a density of 
states. The classical linear approximation for the vibra- 
tional entropy S V ib per particle is then: 
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where N is the number of particles, uj labels the fre- 
quency of the 37V modes, T is the temperature, k B the 
Bolztmann constant. In Eq.([T]) the frequencies u)(T) 
are temperature dependent, both because configurations 
change with temperature, and because non-harmonicities 
are present within a configuration. The latter effect im- 
plies that Eq.([T]) is an approximation, the accuracy of 
which was tested using inelastic neutron scattering com- 
pared with calorimetric measurements in selenium [l9j . 
where it was shown to be accurate throughout the glass 
phase. This leads for the vibrational specific heat C V ib: 
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where (X^))^ = ^^^i^/fiN). Several empirical 
studies have used scattering measurements to estimate 
the vibrational entropy via Eq.([T]) excluding non-linear 
effects by hyper- quenching [l8| or not [III [2(|. Unfor- 
tunately measuring the density of states at various tem- 
perature is difficult, although it is numerically feasible at 
limited viscosities [23|, and has not been done systemat- 
ically for a broad range of fragility. On the other hand, 
it is intuitively clear that the vibrational entropy relates 
to the mean square displacement of the particles (u 2 ) on 
short time scales [l3[ , which indicates the space available 
for particles to fluctuate, and which is well studied ex- 
perimentally. If the particles motions are de-correlated 
from each other, then the volume of phase space in this 
harmonic approximation is: 
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where a labels the three spatial coordinates, <5i?i jQ is the 
displacement of particle i along the direction a, <5Pj jQ , is 
the associated kinetic momentum and m is the particle 
mass. The vibrational entropy is then: 

S mb = k B ln(fi) = h B H{u 2 )) + So + h B ln(T) (3) 

where So is a constant. This expression in general over- 
estimates vibrational entropy, because particle motions 
are correlated: (6Ri, a SRj,p) ^ 0. Eq.Q is therefore 
an upper bound. This expression was discussed in [liij . 
where it was shown to over-estimate at least by a factor 
two the total jump of specific heat of ortho-terphenyl. 
Eq.Q leads for the specific heat in this approximation: 
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The general relation can be obtained by considering 
the linear expression for the mean square displacement: 
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Differentiating Eq.Q with respect to temperature and 
using Eq.@ one gets: 
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Thus Eq.Q, is exact if the effect of temperature is to 
re-normalize all frequencies in the spectrum by the same 
factor, i.e. dhx(ui)/dT is independent of w, leading to a = 
}u) 1 in Eq.©. Nevertheless there are obvious limitations 
in using Eq.Q as it stands to estimate the vibrational 
specific heat in molecular liquids. 

(i) For molecules a fraction of the degrees of freedom 
are very stiff, and do not play a role near the glass transi- 
tion (e.g. the vibrations of strong covalent bonds) . Those 
degrees of freedom do not participate significantly to the 
mean square displacement of the particles. It has been 
proposed to describe molecules as a set of independent 
beads [261] ■ and to estimate the number of beads per 
molecule via the entropy of fusion, compared to the en- 
tropy of fusion in a hard sphere liquid shs — 1.16&B. 
This leads to a rough estimate for the number of beads 
per molecule n^ad = (A_ff /T m )/siis> where T m and AH 
are respectively the fusion temperature and enthalpy. 
Following this line of thought, Eq.Q should be under- 
stood as the entropy per bead, rather than per atom. We 
shall use the notation C V ib to designate specific heat per 
beads, and C V ib for molar quantities. 

(ii) The assumption of overall shift of the typical in- 
teraction stiffness is in general violated. It is well known 
that, at least for fragile liquids, the shape of the spec- 
trum changes with temperature. An excess of soft modes 
with respect to the Debyc model for the density of states, 
the so-called Boson peak (22[, shifts toward higher fre- 
quencies as the system is cooled [H, HU, HH]. Those 
modes presumably contribute significantly to the change 
of mean square displacement and to the vibrational spe- 
cific heat [ill, IHj]. If the increase of (u 2 ) under heating 
is mostly due to the softening of a limited number of 
low-frequency modes, the gain in entropy is diminished 
because the short time scale dynamics becomes more cor- 
related, as can be computed directly from Eq.©. Using 
typical values for the Boson peak will lead to the estimate 
1/5 < a < 1. We shall in fact see below that a must be 
smaller than 1/2 for fragile liquids. Henceforth I shall 
assume a to be constant and independent of the liquid 
fragility, and later come back on those assumptions and 

(4) estimate the effect of the Boson peak on a. 
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To compute the vibrational entropy from Eq.Q, one 
may use the observation that the a-relaxation time scale 
t and (u 2 ) appear to be universally related 0: 



approximation leads to: 
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where (uj> g ) is the mean square displacement at the glass 
transition T g , extracted from the Debye- Waller factor 
measured at the vibrational (picosecond) time scale. In 
some liquids fast relaxation occurs at the nanosecond 
time scale [2{|, which will presumably contribute to the 
total entropy as well. This putative "fast relaxational" 
entropy term is not incorporated in this analysis. Note 
that r and (u 2 ) also appear to be strongly correlated in 
numerical studies of grain boundaries [28|. Eq.® is ob- 
served to hold on a very wide range of liquids fragility, 
the scaling function found is f(x) = fa + fax + fax 2 with 
fa = -0.424, fa = 1.622 and fa = 12.3. Using Eq.© it 
is straightforward to obtain: 
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where Si is a constant. Replacing log(r) by its Volger- 
Fulcher expression log(r) = C + U/(T-T ) in Eq.® one 
obtains an expression for the vibrational entropy which 
diverges logarithmically at To toward — oo (this is obvi- 
ously unphysical and Eq.© or the Volger-Fulcher must 
break down before To). Under the assumption that the 
crystal entropy is sufficiently below the glass entropy near 
T g , an extrapolation of the vibrational entropy and its 
logarithmic divergence must therefore cross the crystal 
entropy at a temperature Tk near To. In this interpre- 
tation of the (Tk, T ) correspondence, there is nothing 
specific about the crystal entropy itself, excepted that 
it is significantly lower than the entropy of the liquid. 
When this is not so, in particular in hard sphere liquids 
where the entropy of the liquid is lower than the entropy 
of the crystal above the melting pressure, Tk — T> (or 
the equivalent relation for packing fractions <f>K = <t>o) 
does not occur- Tk is in fact ill-defined. The relation be- 
tween the apparently distinct temperatures Tk and To is 
generally interpreted as a strong support for the presence 
of a thermodynamic transition where the configurational 
entropy would vanish Q. Nevertheless as soon as a re- 
lation between r and (u 2 ) such as Eq.© holds, those 
quantities are not independent, and an apparent diver- 
gence in r must lead to an apparent fall off in entropy 
at the same temperature, correlating Tk and To- even 
thought no real divergences exist in this interpretation. 

Below T g , one may assume the system to remain in the 
same configuration. Neglecting non-harmonic effects in 
the glass implies (u 2 ) oc T leading to C V i b — 3ks- This 
assumption appears to be a reasonable approximation for 
T/Tg < 0.8 and up to T g in some glasses like polybutadi- 
ene, it nevertheless presumably overestimates the jump 
in specific heat in selenium where non-linearities are still 
significant just below T g Q . Together with Eq. © , this 
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where to = — 91og(r)/91nT|T g is the liquid fragility 
[27| and /'(l) « 26. Ea.lflQ]) thus predicts the jump 
in vibrational specific heat to be strongly correlated to 
fragility. In particular for strong liquids with fragility 
in the twenties, even considering the maximal bound 
a = 1 one finds that AC v a is essentially negligible. Ex- 
perimentally it is found that fragility and molar jump 
of specific heat [16| satisfy (for non-polymeric glasses): 
to = 56T g AC p (Tg)/AH. From Eq.^ we get for the 
slope of the relation between fragility and jump in specific 
heat to = 52AC vlb /(3ak B ) = 52AC vib /(3n bead ak B ) = 
%T m AC vib /AH = ^T g AC mb /AH, where the empir- 
ical rule T m = 3/2T g was used. Thus if the jump 
in specific heat was fully vibrational one would have 
a = 0.53 ~ 1/2. For fragile liquids a cannot be sig- 
nificantly larger than that value, otherwise the jump of 
vibrational specific heat would be larger than the total 
jump, which must also contain the positive contribution 
of the configurational entropy. 

I now come back to the value and dependence of the 
parameter a, which contains both the overall stiffening 
of the structure, and the softening of a limited num- 
ber of degrees of freedom, the Boson peak. To esti- 
mate a in the case where only the Boson peak affects 
the mean square displacement, I consider the follow- 
ing approximation: the density of states is assumed to 
follow the Debye law D(uS) cx w 2 /w|j for u < lobp, 
where ujd is the Debye frequency. D(uj) remains of the 
same order D(uS) cx 1/ujd for uibp < w < Wd. This 
corresponds to a Boson peak (maximum in D(uj)/lu 2 ) 
at uibp- If wbp(T) is significantly smaller than wj, 
one finds for such a spectrum (l/w 2 ) w « 1/(udUbp). 
Furthermore I assume that the modes near the Boson 
peak are those sensitive to the temperature, leading 
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tice dln(u}Bp)/dln(T) < (the system softens with tem- 
perature), and one gets a > ujbp/ojd- Typically values 
of the Boson peak in the glass are ujbp ~ &d/5 
leading to a > 1/5. From this primitive estimate of the 
lowest bound of a, one also gets that a increases as the 
system is cooled, as does the Boson peak frequency, and 
is presumably smaller in fragile liquids where the shift 
of the Boson peak with temperature is pronounced. If a 
increases under cooling, the apparent divergence of the 
extrapolated vibrational entropy will only be enhanced, 
and the interpretation of the To = Tk correspondence 
remains safe. On the other hand, the possible system- 
atic decay of a with fragility may diminish the predicted 
correlation between the vibrational part of the jump in 
specific heat and fragility. Taking our lowest bound of 
a R3 0.2 one still finds the vibrational contribution to the 
jump of specific heat to be significant (about 1/3) for the 
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most fragile liquids. This is consistent with the observa- 
tions that the vibrational fraction of the jump of specific 
heat increases with fragility and represent more than half 
of the jump for fragile liquids (2lj . 

I have argued that as soon as a relation between relax- 
ation time and mean square displacement such as Eq.([8j 
exists, dynamics and vibrational entropy must be cor- 
related, leading to an alternative explanation for the 
Tq = Tk correspondence. The amplitude of the corre- 
lations between fragility and jump in specific heat im- 
plies that the mechanism at play must include the soft- 
ening under heating of a limited number of soft degrees 
of freedom, and is not simply a change of elastic mod- 
uli. This likely corresponds to the well-known softening 
of the Boson peak, which has been shown experimentally 
to affect significantly the thermodynamics of some liquids 
[l8j . Assuming that this is the case I find a reasonable es- 
timate for the magnitude of the vibrational specific heat. 

It has been proposed since Adam-Gibbs [7H9j that the 
relation between excess entropy and dynamics observed 
in liquids indicates that as less configurations are visited, 
more collective rearrangements are required to relax the 
liquid, causing a rapid slowing down of the dynamics. 
However, there appears to be a viable alternative possi- 
bility consistent with this observation: when a liquid is 
cooled favored structures appear. This lowers the con- 
figurational entropy but more importantly stiffens the 



liquid which causes activation barriers to increase and 
vibrations to decrease. In this "elastic" view of the glass 
transition the vibrational entropy most directly reflects 
the liquid dynamics, in consistence with the present anal- 
ysis. Empirical evidences for such a stiffening range from 
vibrational spectra [H, H3, l25l- mean square displace- 
ments and elastic moduli[13. l3l| measurements. In 
all cases strong correlations with the dynamics are found. 
Consistent with this approach, liquids at constant volume 
which are known to stiffen significantly less fig] . [231 than 
at constant pressure are also much less fragile [32| in the 
cases known. Theoretically, the view that the stiffening 
of soft modes control the dynamics is supported among 
other works by the mode-coupling theory at moderate 
33L 33 | , by geometric arguments on the nature 



viscosities 

of the Boson peak and pre-vitrification in hard sphere 
liquids [HI and by the numerical observations that ac- 
tivated events required to flow occur mostly along soft 
modes belonging to the foot of the Peak [36j, [37j ■ Theo- 
retical efforts required to assess this view further should 
focus in particular on the still obscure relation between 
excess soft modes and viscosity in highly viscous liquids, 
and on the questions of what in its microscopic structure 
determine if a liquid stiffens or not under cooling. 

It is a pleasure to thank C.A. Angell, M. Cates, G. 
Tarjus and F. Zamponi for comments on the manuscript. 
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